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Rapidly rotating Bose-Einstein condensates confined in anharmonic traps can exhibit a rich variety 
of vortex phases, including a vortex lattice, a vortex lattice with a hole, and a giant vortex. Using 
an augmented Thomas-Fermi variational approach to determine the ground state of the condensate 
in the rotating frame - valid for sufficiently strongly interacting condensates - we determine the 
transitions between these three phases for a quadratic-plus-quartic confining potential. Combin- 
ing the present results with previous numerical simulations of small rotating condensates in such 
anharmonic potentials, we delineate the general structure of the zero temperature phase diagram. 

PACS numbers; 03.75.Fi, 67.40.Db, 67.40.Vs, 05.30.Jp 



I. INTRODUCTION 

Gaseous Bose-Einstein condensates of alkali-metal 
atoms are fertile systems in which to study the behavior 
of superfluids under rotation. For rotation frequencies 
above Sid , the critical frequency for creating a single vor- 
tex, the properties of these systems resemble those of liq- 
uid helium under typical experimental conditions. States 
with a small number of vortices have been observed in 
the experiments of Refs. As increases, more and 

more vortices form, and indeed it has become possible 
to create lattices [^,|| of more than 100 vortices. The 
physics of rapidly rotating gases is different in harmonic 
and anharmonic traps. In a harmonic trap, the rotation 
rate is limited by the radial trap frequency w, and sets 
the scale of fici ■ As the angular frequency of rotation of 
the gas approaches the transverse trapping frequency, the 
centrifugal force becomes equal to the restoring force ex- 
erted by the trap and the atoms are no longer contained. 
On the other hand, the regime Vl> uj becomes accessible 
in anharmonic trapping potentials that rise more steeply 
than quadraticaUy. 
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FIG. 1. Schematic phase diagram of the ground state of a 
rapidly rotating Bose condensed gas at zero temperature in 
the Q.-oas plane, showing the regions of multiply quantized 
vortices, of singly quantized vortices forming an array which 
at large number becomes a triangular lattice, and of an array 
of singly quantized vortices with a hole in the center. The 
figure does not sketch out the details of the phase diagram at 
small S7 below the critical rotation, fici, for onset of vorticity. 



Condensates rotating in anharmonic trapping poten- 
tials are expectd to have a rich structure of vortex states, 
as described in Refs. [HJIJ^,^, reviewed below. These 
states include an array of singly quantized vortices (which 
becomes a triangular lattice when the number of vortices 
is large) ; a vortex array with a finite radius "hole" in the 
center - opened by the centrifugal force ~ in which the 
condensate density is zero but the vorticity non-zero; and 
giant (multiply quantized) vortices. Our aim in this pa- 
per is first to study, via a variational approach, the states 
in an anharmonic quadratic-plus-quartic radial potential 
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where p is the (cylindrical) radial coordinate, m is 
the atomic mass, the oscillator length 

{Ti = 1), and A is a dimensionless measure of the an- 
harmonicity. We consider a container infinitely long and 
cylindrically symmetric and rotating about the z axis, 
with a density per unit length a — N / Z in the axial (z) 
direction; the particles are assumed to interact via s-wave 
scattering with scattering length as > 0. The variational 
approach is adopted from Ref . jsj , with the density profile 
of the gas given by an augmented Thomas-Fermi approx- 
imation which explicitly takes into account the important 
kinetic energy of flow about the individual vortices and 
the modification of the interaction energy by the den- 
sity variations due to the vortices. We then combine the 
present calculation with the results of the earlier studies 
to delineate the phase diagram of a rotating Bose conden- 
sate in an anharmonic trap, when the total vorticity is 
3> 1 (i.e., well above the critical rotation frequency f2ci)- 
As we shall see, the phase diagram has the schematic 
structure in the fl-aug plane shown in Fig. 1. At very 
small interaction strength or density multiply quantized 
vortices are always preferable, as found in f^, while at 
larger interaction strength at small f2 the system forms 
an array of singly quantized vortices. At larger the 
array, as shown in |Q, develops a hole in the center. The 
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phase diagram contains a triple point between the three 
regions. 

Lundh by numerical simulation of the Gross- 
Pitaevskii equation in two dimensions, determined the 
ground state vortex structures in a rapidly rotating con- 
densate in radial trapping potentials ^ p" (n > 2), as 
well as in the potential (|l|) for small total vorticity, ^ 10. 
His calculations show that for small dimensionless inter- 
action strength, aag, the system prefers multiply quan- 
tized vortices, as illustrated in Fig. 1, while above a crit- 
ical (A-dependent) interaction strength, {aas)c ~ 1, the 
ground state consists of an array of singly quantized vor- 
tices. The fact that a giant vortex state is favored at 
small fTfls can be understood from the fact that in the 
limit of a non-interacting gas trapped in a potential that 
rises faster than quadratically, the dependence of the sin- 
gle particle energies, eji'), on angular momentum hv rises 
faster than linear P,p^. Thus, in this limit, in the ground 
state of a fixed number of bosons in the frame rotating 
at given fl all the particles occupy the single angular mo- 
mentum state that minimizes e(i/) — flv. 

Condensates in the potential (Q) were also studied via 
numerical simulation by Kasamatsu, Tsubota, and Ueda 
[0, at fixed parameters A = 1/2 (corresponding to their 
k = 1) and aug — 250/87r. At low H. the system con- 
tains an array of singly quantized vortices. [The value of 
aUg in this calculation well exceeds the critical coupling 
{aas)c found by Lundh, and is above the triple point in 
Fig. 1.] However, as they increase the rotation frequency 
to 11 ~ 2.14a;, the vortex array begins to develop, as a 
consequence of the centrifugal force, a hole in the center. 
At still higher rotational velocity, ^l ^ 3.2 — 3.5, the sys- 
tem undergoes a transition to a single multiply quantized 
giant vortex. 

Similar behavior was found in Ref. |^ for a condensate 
contained in a cylindrical hard walled bucket, by con- 
structing an augmented Thomas-Fermi solution of the 
Gross-Pitaevskii equation. At low rotation speed the vor- 
tices form a uniform lattice. However at a critical rota- 
tion frequency, l/2m^o^ ~ where ^0 = (87rnas)~^^^ 
is the coherence or healing length, and n is the average 
density - the lattice begins to develop a hole in the center. 
Finally at a second critical frequency, JIg ~ l/2w^0) the 
system makes a transition to a giant vortex state. Fet- 
ter 1^, in considering condensates in the potential (||), 
finds a critical rotation speed at which the density in the 
center of the trap goes to zero; however, the onset of 
the holethe is obscured by the particular approximation 
employed following factorization of the condensate wave 
function. 

While detailed comparison of our calculation with the 
numerical simulations of Refs. |^ and |^ is hindered by 
the fact that the number of vortices in the simulations is 
not large enough to form a uniform lattice, as we assume 
(see, e.g.. Fig. le of Ref. |Q, and Fig. 6 of Ref. our 
results are qualitatively consistent with and provide a 
simple theoretical understanding of the numerical results. 
At large fl/co and aug where our approach is valid, we 



see that with increasing the system indeed undergoes 
a transition from a vortex lattice to a vortex lattice with 
a hole around the axis of rotation. Furthermore, as f2 
increases the system eventually undergoes a transition 
to a single giant vortex, for aas up to at least 350, the 
largest value for which we studied this transition. 

We first describe the three states of the rapidly rotating 
gas - a vortex lattice, both with and without a hole along 
the axis, in Sec. IIA, and a giant vortex, in Sec. IIB - and 
calculate the augmented Thomas-Fermi solutions for the 
lowest energy in the rotating frame at given Q. In Sec. Ill 
we determine the transitions between different phases for 
the particular anharmonicity A = 1/2, and the less steep 
A = 1/8. In Sec. IV we examine the limits of validity of 
our approach, and finally summarize our results in Sec. V. 



II. POSSIBLE VORTEX STATES OF A RAPIDLY 
ROTATING GAS 

A. Vortex lattice 

We first consider a uniform lattice of singly quantized 
vortices. As shown in Ref. ||], the energy of the system 
in the rotating frame, E[^, in the Thomas- Fermi limit is 



EL- 



1, 



q2 

- ) mp^ + (a + i)rj„ - VL 



(2) 



where is the rotational frequency of the trap, r2„ 
is the rotational frequency of the vortex lattice, and 
= ATrh^as/m. We assume an effective core radius 
^ (different in general from the Gross-Pitaevskii healing 
length ^o)- The function a(C) describes the kinetic en- 
ergy of rotation of the vortices, and &(C) the effects of 
the vortex cores on the interaction energy. For a simple 
model in which the order parameter ramps up linearly 
over a distance ^ from the cores of the vortices, and then 
flattens out, a and b are given by 
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where C = and i = l^/mTl^ is the radius of 

the cylindrical Wigner-Seitz cell corresponding to a given 
vortex 

Minimization of with respect to ni^{p) yields the 
Thomas-Fermi profile, 
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with p the chemical potential. Since the density of the 
gas at the center of the cloud, p = 0, cannot be negative, 
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the system begins to develop a hole when tt-l (p = 0) falls 
to zero, or from Eq. /i — (a + l)ri„ + = 0. For 
smaller fl, this quantity is positive and the inner radius 
of the cloud, vanishes, while for larger ft, Ri ^ 0. 
We note that when a hole develops in a uniform vortex 
lattice, the hole must, for the lattice to remain uniform, 
contain the same vorticity, 7ri?f — mRifly (where Uy is 
the number of vortices per unit area) , as would be present 
were the hole filled with condensate. Equation (^) implies 
that the inner and outer radii of the cloud, Ri and R2, 
respectively, where the density vanishes, are given by the 
solutions of 
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where here lengths with a tilde are measured in units of 
d, e.g., Ri = Ri/d, i = \,2, and energies or frequencies 
with a tilde are measured in units of uj. 



Integrating the density over the volume, J riL (P. 
we find, 

^aa,b = fi{Rl-Rl)--^(Ri-R\) 

-^{Rl - Rl) + [Rl - Rl)[^ - (a + mv] 
-\{R'2-R\Wl-'i^^v). 

Minimization of E'^ with respect to Vt^ yields 
Vty ^Q- N{a+l)/I, 
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where I = J (Pr mp^ni^{p) is the moment of inertia. 
Equation , with Eq. (^) , implies that 
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For a given value of the dimensionless parameters eras, 
A, and (, we solve Eqs. (||), (|), and (||) to find i?2, 
/i, and fly-, then from Eqs. (|) and (Q) we construct the 
energy £"^($1,0 in the rotating frame, as a function of 
fl and Cj and finally vary to minimize Ei. The total 
angular momentum of the gas is given by pi[ 
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A useful check on the numerical results is the relation 
dEl/dn = -L„ implied by Eq. (0). 



B. Giant vortex 

We next consider the giant vortex state, whose vortic- 
ity is localized along the rotation axis, and order param- 
eter is 
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where nc^p) is the density profile of the cloud, p and (f) 
are cylindrical polar coordinates, and v is the winding 
number. The angular momentum per particle along the 
z axis is Lz/N = v, while the energy of the system in the 
rotating frame in the Thomas-Fermi limit is 
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Setting the variational derivative 5Eq/ 5nQ(p) equal to 
the chemical potential, /x, we find the density profile of 
the cloud: 
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where p' = p + vU,. The inner and outer radii of the 
cloud, Rl and i?2, solve 
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In addition, the constraint J uq d^r — N implies that, 
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Minimization of Eq with respect to v yields 
no 
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and thus 



2~p'\n^ 

Rl 



[Ri-Ri) 



V 

~2 

_ A 
' 4 



1 1 

^2 ^1 
t^2 ~ -f^l 



(16) 



For given values of the dimensionless parameters aa^ and 
A, we solve Eqs. (p^, (p^, and (|l6|), and express i?2, 
/i' and O in terms of v. Wc then numerically solve for 
E'q from Eqs. ( |ll| ) and (|l^) as a function of v, and thus 
implicitly as a function of $7. The relation dE'^jd^ = 
— Lj, implied by Eq. (pl|), again provides a check on the 
numerical results. 
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III. RESULTS 

We have limited the calculations of the vortex phases 
here to the cases of intermediate anharmonicity, A = 1/2, 
and weak anharmonicity, A = 1/8, which well illustrate 
the richness of the phase diagram. Figure 2 shows the 
phase diagram for A = 1/2. Relatively small ft and rel- 
atively large aug favor the formation of a vortex lattice. 
In the opposite limit, large f2 and small aas favor the 
formation of a giant vortex. Finally at intermediate f2 
(for aas in the range considered) the system prefers a 
vortex lattice with a hole. Quite generally, the transi- 
tion from a vortex lattice to a vortex lattice with a hole 
is continuous, second order, while the transition from a 
vortex lattice with or without a hole to a giant vortex is 
discontinuous, first order. 
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FIG. 2. The phase diagram of the rapidly rotating Bose 
condensed gas at zero temperature in the fl-aas plane, for 
anharmonicity A — 1/2. The labels VL, VLH and GV denote 
the regions where the favored state is a vortex lattice, a vortex 
lattice with a hole, and a giant vortex, respectively. The filled 
circles at aa, > 200 indicate the interaction strengths explic- 
itly calculated here, while the two open circles at aas ~ 10 
show the calculations of Ref . [7] . The lines are a guide to the 
eye. 



in the calculation) is between 2 and 3. For such low 
values of v our approach is clearly not valid, and the 
discrepancy should not be taken seriously. 





-o 




-6 


LU 


-7 


s" 


-8 






LU 


-9 




-10 




2.34 2.37 



2.43 



2.4 

FIG. 3. The energies in the rotating frame of the giant 
vortex, E'q (dashed curve) and of the vortex lattice, (solid 
curve) as a function of Q, for anharmonicity A — 1/8, and 
dimensionless interaction strength aas ~ 90. 
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FIG. 4. The angular momentum per particle, Lz/N, versus 
the angular frequency Q, for A = 1/8 and aas ~ 90. 



For (Tfls = 200, the vortex lattice develops a hole for 
U,/lo K 3.90, and then makes a transition to a giant vortex 
for Vi/ix) ~ 5.48. At this point the giant vortex has a 
vorticity v ~ 168. For aag = 550, corresponding to the 
parameters of the multivortex MIT experiment Ref. , 
we find that the vortex lattice develops a hole for Vl/uo k, 
4.50, and that even up to the highest value oi^l/uj k, 8.4 
that we examined, this solution has lower energy than 
the giant vortex. For such high values of eras and we 
observe that the two solutions (vortex lattice with a hole 
and giant vortex) look macroscopically very similar. 

For comparison, we have also calculated the phases for 
interaction strength aas = 250/87r « 10 considered in 
Ref. , although it is well below the limits of validity 
of the present approach. Then the transition to a giant 
vortex occurs for ^/uJ ~ 1.5, compared with the value ^ 
3.2 found in Ref. 0; however v (taken to be continuous 



The particular case of weak anharmonicity case A — 
1/8, with aas — 90 also shows the three phases: for 
Q,/u} K, 2.30 the vortex lattice develops a hole in the 
middle, while for Vt/u) 2.40 it makes a transition to 
a giant vortex. Figure 3 shows the gentle crossing of 
i?Q(il) (solid curve), and E'^{^) (dashed curve), while 
Fig. 4 shows the angular momentum per particle L^/N 
(of the lowest state) , the negative of the slope of the lower 
of the two curves of Fig. 3, as a function of Q.. Since the 
curves cross at a finite but small angle, L-^/N is discon- 
tinuous at the point of the transition. Figure 5 shows Ri 
and i?2 as a function of Q. for the three solutions. 
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FIG. 5. The inner and outer radii, Ri and R2, of the cloud 
vs. the angular frequency il, for A = 1/8, and interaction 
strength aas = 90. The solid curves are for the vortex lattice, 
and the dashed ones for the giant vortex. As the lowest curve 
illustates, the vortex lattice develops a hole for Q/uj ~ 2.30. 
The transition to the giant vortex takes place at Q/uj « 2.40. 



IV. LIMITS OF VALIDITY OF THE PRESENT 
APPROACH 

The present approach makes two crucial assumptions, 
first that the profile of the gas can be described by the 
Thomas-Fermi approximation, with the vortex correc- 
tions to the kinetic and interaction energies included, 
and second that the array of singly quantized vortices 
is a uniform lattice. The validity of the Thomas-Fermi 
approach requires that the kinetic energy per particle, 
~ l/2TO(i?2 — ^1)^, associated with localizing the atoms 
within the radii of the profile of the cloud Ri and R2, 
be much less than nllo = 1/2to^q, where n ~ N/tt{R2 — 
Rl)Z is the average density of the cloud, and ^0 is the co- 
herence length. Thus validity of Thomas- Fermi requires 
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or equivalently, <^ R2 — Ri- For Ri = 0, the condition 
( |l7| ) is simply Saag » 1. 

In addition to condition (p7|), the basic features of the 
lowest state have to be captured by the trial functions 
being considered. For small values of aus and Q/uj, the 
numerical solutions have too small a vorticity to 
arrange in a triangular lattice. Rather, in Ref. |Q at 
f2/u; = 2.5 for example, the vortex hole, containing 4 
units of vorticity, is surrounded by a ring of 10 vortices, 
instead of a developed lattice. A reasonable estimate 
for the lowest value of at which the calculation is 
valid is ~ 50, which typically (e.g., at A = 1/2) requires 

n/io > 3.5. 



V. SUMMARY 

In the present study we have examined, by means of 
a variational method, the lowest state of a Bose-Einstein 
condensate rotating in a quadratic-plus-quartic confining 
potential, in the limit of fast rotation. We find that as 
the frequency of rotation, the density of the gas, and the 
strength of the quartic term in the trapping potential 
are varied, the system undergoes transitions among three 
phases: a vortex lattice with and without a hole, and a 
giant vortex. The phase diagram, whose structure we 
have mapped for a specific choice of the coefficient of 
the anharmonic term, is quite rich. This structure is not 
limited to the specific form of the anharmonic potential 
considered here. Our approach is reliable in the limit 
of large densities and high angular frequencies. Many 
problems remain, including determining the details of the 
critical point in Fig. 1, the variation of the phase diagram 
as the coefficient of the quartic term in the confining 
potential is varied, and finally, the structure of the phase 
diagram at finite temperatures. 
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